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ABSTRACT 

We study a model for two-dimensional skyrmions on a sphere of radius L. 
Such model simulates a skyrmion lattice of density W/{2'n'L'^), where W is 
the skyrmion winding number. We show that, to a very good approximation, 
physical results depend only on the product aL^, where a is the strength of 
potential term. In the range aL^ < 3 the order parameter vanishes, there is a 
uniform distribution of the density over the whole surface and the energy of the 
W = 2 sector lies above twice the energy of the W = 1 sector. If aL^ > 6 the 
order parameter approaches unity and the density concentrates near one of the 
poles. Moreover the disoliton is always bound. We also present a variational 
solution to the field equations for which the pure aL^— dependence is exact. 
Finally, some consequences of our results for the Quantum Hall Effect are 
discussed. 
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Although skyrmions were originally introduced to describe baryons in the context of 
non-linear meson theories [|1| , they are playing an increasing role in other areas of physics 
as well. In the present paper we consider the case of "baby skyrmions", i.e. skyrmions in 
two spatial dimensions. For example, 2D skyrmions are used to model the bubbles that 
appear in certain condense matter systems in the presence of an external magnetic field ; 
they could provide a mechanism associated with the disappearance of antiferromagnetism 
and the onset of HTc superconductivity; etc.. Recently it has been suggested that they 
may be responsible for many features of the Quantum Hall Effect (QHE). Baby skyrmions 
have been also studied in the context of strong interactions as a toy-model in order to 
understand the more complicated dynamics of usual skyrmions which live, of course, in 3 
spatial dimensions 0. 

For many of these applications it is of great interest to consider not only the proper- 
ties of few baby skyrmion systems but also the behavior of such properties as the density 
of skyrmions increases. So far, this problem has been studied by putting baby skyrmions 
on a lattice, which is, of course, a rather cumbersome numerical task (see e.g. Ref. [Q). 
Analogous studies have been done for 3D skyrmions to understand baryon properties in 
dense matter and the appearance of phase transitions such as chiral restoration, decon- 
finement, etc [^. However, in the 3D case it was realized that many of the qualitative 
results of the lattice calculations (and even some quantitative ones) could be obtained in 
a much easier way by studying the behavior of few skyrmions on a hypersphere 0. This 
approach was introduced in Ref. and proposes to replace the periodic array in flat 
space by a one (or few) skyrmion system on the compact manifold S^(L) of radius 
L. The finite baryon number on S^(L) corresponds to a finite baryon density, so that 
one obtains a model for skyrmion matter which is much simpler to study than any lattice 
model. The density of this matter can be increased by decreasing L. The aim of the 
present paper is to extend these ideas to the case of 2D skyrmions on the sphere. 

Baby skyrmions are obtained as non-trivial solutions of the well-known non-linear 
0(3) model. This model consists of three real scalar fields (pa (o = 1,2,3) subject to the 
constraint (p ■ (p = 1- The equations of motion admit static solutions with finite energy 
which represent a mapping of H^pat iiito Sf^^. They are characterized by the winding 
number W and the density p 

W = ^Jd^rp ; p = e,, 0-(9>x^-0) (1) 

where i,j stand for the spatial coordinates. We start from the energy functional 

E = ^(2) + ^(4) + 

E^^^ = ]^j(frdS-d'$ E^'^ = ^d'rpp- E^^^ = j dh {n^ ~ (2) 
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where 77,3 is a unit vector in the third direction in internal space and a is a parameter 
that is assumed positive. This type of energy functional is obtained as the static limit of 
either the relativistic 0(3) model or the non-relativistic Chern-Simons theory. 

In 2D the quadratic term in the field derivatives gives a scale independent contri- 
bution to the soliton mass, due to the well known compensation between the derivatives 
and the integration pj. Thus, the stabilization procedure for baby skyrmions is usually 
based on a competition between a quartic term and a linear potential term. This last 
(attractive) interaction tends to concentrate the skyrmion at 03 = 1, yields a contribution 
to the total energy proportional to the square of the soliton size A and destroys the orig- 
inal 0(3) symmetry. In condense matter physics this potential is usually related to the 
magnetic interactions (e.g. Zeeman term in QHE). In the framework of hadron physics 
it would correspond to the explicit chiral symmetry breaking driven by the pion mass 
term and therefore considered to be only a small perturbation, and often neglected. The 
collapse of the field is prevented by means of the quartic term which yields a contribution 
proportional to In QHE this contact interaction would be a "poor man" mocking 

of the Coulomb repulsion, while in hadron physics it can be understood as a large mass 
limit of the p-meson exchange. It should be mentioned that, as in the 3D case, the quartic 
term can be replaced by higher order derivative termsQ. From the experience in 3D we do 
not expect major changes in the behavior of most of the quantities to be discussed below. 

The simple relation between the coefficients of the quadratic and quartic terms in 
(0) is obtained through a convenient choice of the unit of length. Therefore, in flat space 
there is only one parameter in the model, namely the strength a of the potential term, 
aside from an overall multiplying factor which is already omitted in (^. 

It is interesting to note that using the inequality 

< \d\ ^(9,0±e,,0xa^0)%^Qe,,9>x9^0±v^(n3-0))' (3) 

one can find the Bogomol'nyi bound for the baby skyrmion energy. One obtains 

E > 4:7ik (1 + v^) (4) 



which improves over the usual bound through the contribution proportional to y/a. 

As mentioned above, our aim is to extend the model above by going from R^^ to 
Sgp(L) where L is the radius of the two-sphere. A convenient choice of coordinates on 
S^p(L) is given by the conventional polar coordinates 9,ip {0 < 9 < n and < (p < 2n) 

x = Lsm6cos!f; y = Lsm6sm(p (5) 



^In fact, using particular combinations of derivative terms it is possible even to stabilize the soliton 
in the absence of the potential term. Although interesting, we do not consider this possibility here. For 
details, see Ref.lH. 
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The jacobian of the transformation and the metric associated with the polar coordinates 
are 



(6) 



In order to obtain exphcit static solutions in the winding number W = k sector we 
introduce the hedgehog parameterization 



sin / cos kip 



sin / sin k(p 



cos/ 



(7) 



where the "radial" profile / = f{6) is subject to the boundary conditions /(O) = tt and 
/(vr) = 0. 

Using (0) and we may express E^, the soliton energy divided by Arck, in terms 
of the profile / as 



= — i de sin e 

Ak 



'sin/\ ^ _P / sin/ 
sin 6 I L"^ \ sin 6 



+ 2aL2(l - cos/) 



while the winding number density results 



Pk 



— f ^^^^ 
sin 9 



(9) 



In Eqs.(H,|^) a prime denotes differentiation with respect to 6. 

Minimizing we obtain the eq. of motion for /. It results 



1 + 



sin/' 



LP' \ sin 9 

+ 



J, _ sin 2/ ^ fc^^; Sin/ /^, cos/ _ sin/ ' 
sin 26* sin^ V cos 6' sin^ 



cote - aL^sin/ = (10) 



Fig. 1(a) displays the profile function f{9) in the k = 1 sector for different values of 
L, and a = 0.1. Fig. 1(b) shows the same magnitudes but for different values of a, and 
L = 2. For very dense systems or for very small values of a the identity map f = n — 9 
constitutes a valid (although not exact) solution for /. As shown in Figs. 1(c) and 1(d), 
it corresponds to a uniform density pi over the surface of the sphere. As L or a increases, 
both the profile and the density become more concentrated near ^ 0. 

The fact that similar results are obtained by varying either L or a raises the question 
whether the baby skyrmion on the sphere is a single-parameter model. The origin of such 
uniqueness may be traced back to the fairly good independence of the quadratic term E^^ 
respective to both L and a . This can be clearly seen in Figs. 2 (a) and 2(b). There it is 
shown the behaviour of the different contributions to the soliton energy as function of L 
and a, respectively. As a consequence, the physics of the model is determined only by the 
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1/2 

(12) 



interplay between the quartic and the potential term, and thus depends basically on the 
product aL"^. The fact that some quantities (e.g. the mean square radius in Fig.2) behave 
differently as function of L or as function of a is simply due to explicit dependencies on 
L. 

It is clear from Fig.2 that we may distinguish two regions which display different 
behaviors of the physical magnitudes. For large values of aL'^ the order parameter defined 
as the mean value of 03 

^■, = ^1 desine (I), (11) 

tends to 1; the quartic and potential contributions E^^^ < E^^^ become rather similar; the 
energy of the k = 2 sector is smaller than (twice) the one of the k = 1 sector; the mean 
square radius 

ffc = 2tcL'^ J d9 sin 9 pk (sin i 

decreases, if plotted in units of L. On the contrary, for small values of aL^ the order 
parameter vanishes; E'*^^^ >> E'^^^; the energy of the k = 2 sector is larger than (twice) the 
one of the k = 1 sector; the mean square radius remains approximately constant. Between 
the two extreme behaviors there is a broad transition region corresponding approximately 
to the values 3 < aL'^ < 6. This differs from the situation in the absence of the quadratic 
term E^^^[^. 

The behaviour of the physical quantities as function of L is, in general, similar to 
the one obtained in 3D 0| . Due to the presence of the potential term, however, the sharp 
transition found in 3D (in the chiral limit) turns out to be a broad one here. Such term 
is also responsible for the absence of "swelling" of the mean square radius in our model. 
This seems to be the only qualitative feature of 3D skyrmions which is not shared by baby 
skyrmions. 

In the context of QHE there is an undergoing discussion about the stability of the 
disoliton at low densities, with respect to the decay into two k = 1 solitons. In Ref. ||TT| , 
it was found that the disoliton is stable independently of the strength of the potential 
term. This seems to be in contradiction with Ref. ||T^, where it is shown that there 



is a critical strength above which the disolitons become unstable. Our results indicate 
that at zero density the disoliton is always stable. They also confirm the conjeture made 
in Ref. ||Tl| that, as density increases (and therefore L decreases), there might be a 
transition above which disolitons become unstable. The fact that at zero density (flat 
space) the disoliton is stable is in agreement with the result obtained in 3D. There the 
disohton is supposed to represent the deuteron. As well-known, at classical level, the 
lowest diskyrmion (the doughnut) is rather strongly bound |T^. Recently |]T4| it has 
been shown that quantum effects tend to decrease such binding. It is clear that in order 
to make more solid statements such quantal effects should be introduced in our model 
through a proper quantization method (see, for instance, Ref. [|^). Another point 
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to take into account when applying our model to QHE is the long range nature the 
Coulomb interaction. As mentioned above a very crude approximation has been used in 
the present work. Namely, we use a local density- density interaction as a four- derivative 
term. The consequences of replacing such term by a long range interaction is currently 
under investigation. 

We finish by presenting an approximate variational solution for the baby skyrmions 
on the sphere. For this purpose one can use the solutions of the equation of motion ( |lO|) 
in the absence of quartic and potential terms. They can be written as 

k 

(13) 



tan — 
2 



7] 9 

- cot - 
2 2 



where rj is an arbitrary constant. These solutions reduce to the known expression in fiat 
space [H and saturate the corresponding Bolomol'nyi bound (i.e. Eq.(^) with a = 0) for 
any value of rj. Therefore when used as variational solution of the full model the value of 
7] turns out be only a function of the combination aL^. For example, minimization of the 
soliton energy in the k = 1 sector leads to the following equation for t] 

a.L' = - ^'''"^ + ^ (14) 

All the features that are present in the exact solution of the hedgehog ansatz are qual- 
itatively reproduced, the quantitative agreement lying within few percents. A detailed 
comparison will be given elsewhere. 

To conclude, we have studied the consequences of locating one or two baby skyrmions 
on the surface of a sphere of radius L. Thus the model covers the whole range between 
the two limits of low densities (flat space) and of high density (small L). The soliton is 
stabilized through the presence of a (repulsive) quartic term in the field derivatives and of 
an (attractive) potential term proportional to the strength a of an external field, in addi- 
tion to the usual quadratic term. A value of the Bogomol'nyi bound which is higher than 
the one currently used in the literature on baby skyrmions has been obtained. Within the 
hedgehog approximation, the model was solved exactly. Although in principle the model 
has two independent parameters, we show that the product aL^ determines the physically 
meaningful properties of the system. The solution displays two characteristic regimes: the 
order parameter vanishes and there is a uniform distribution of the density over the whole 
surface for the range aL'^ < 3 ; on the contrary, the order parameter approaches unity 

and the density concentrates near ^ ^ if aL^ > 6. In the first region the energy of the 
k = 2 sector lies above twice the energy of the k = 1 sector. The opposite is true for the 
soliton region, in which the disoliton is always bound. Finally, a variational solution that 
realizes exactly the pure aL^— dependence is given. We believe that some of the results 
obtained in the present work can be directly extended to interesting issues such as phase 
transitions in Quantum Hall Effect. There, the skyrmion density would be related with 
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the departure of the filhng factor from one. However, some improvements, such as e.g. 
inclusion of long range interactions, have to be made before drawing stronger conclusions. 
We hope to be able to report on these matters in forthcoming publications. 
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Figure 1: Soliton profiles (a) and (h) and winding number densities (c) and (d). The 
densities ^ are plotted in units of '^t^ / L?' . 
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Figure 2: Various properties of the solution. Figs, (a) and (b) contain the different 
contributions ^ to the soliton energy for winding number k=l (in units of Air). Figs, 
(c) and (d) include the total soliton energies the order parameter and the mean 



square radius (|7^j. 
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